Current— phase relation in Josephson junction coupled with a magnetic dot 
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The current-phase relation for a short Josephson junction placed in the nonuniform field of a 
small ferromagnetic particle is studied. The effect of the particle produced on the junction appears 
to be strong due to the formation of the pair of oppositely directed Abrikosov vortices which pierce 
the thin film superconducting electrode and cause a small-scale inhomogeneity of Josephson phase 
difference. The induced phase difference inhomogeneity is shown to result in the nonzero fixed phase 
drop tpo across the junction. The equilibrium value ipo corresponding to the ground state of the 
junction depends the configuration of the vortex-antivortex pair. The possibility to tune the ground 
state phase difference ipo is discussed. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Usually the current-phase relation (CPR) in Joseph- 
son junction close to the critical temperature is sinu- 
soidal Is{v) = -fcSin(p, and the dependence of the 
free energy Ej = {hlc/2e){l — coaip) assumes posi- 
tive values of the critical current Ic > (see Ref. [l|). 
So, in the absence of a supercurrent, Is = 0, the 
phase drop across the conventional junction equals zero^. 
But under certain conditions one can fabricate so-called 
Josephson tt junction^ which has an energy minimum 
at = TT, i.e., it provides a phase shift of tt in the 
ground state (Refs. lUU). The CPR of tt junctions 
reads Isi^fi) = Ic sin((^ -|- tt) and can be formally de- 
scribed by the negative value of the critical current Ic- 
The TT states have been observed in Josephson junctions 
consisting of two d— wave superconductors^ in super- 
conductor/ferromagnet/superconductor (SFS) junctions 
utilizing ferromagnetic barriersSili^, and also in multi- 
terminal superconductor /normal-metal / superconductor 
(SNS) Josephson junctions^. Such tt junctions are sup- 
posed to open up new opportunities for designing Joseph- 
son effect-based devicesiSiiii'i^. 

Recently the investigations of Josephson ip junctions 
which provide the realization of an unusual current-phase 
relation 



Is{ip) = Ic sin((/? -I- ipo) 



(1) 



have been attracting a lot of attentior^i^iiiiii^. The min- 
imum of the Josephson energy of (p junctions Ej — 
(fi/c/2e)[l — cos(</3-|-t^o)] corresponds to the nonzero value 
of the phase difference (p — —ipo such as < (/Jq < tt- The 
realization of the tp junction is possible in the case of pe- 
riodic structures composed of alternating and tt mini- 



junctioni 
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, or in the case of SNS structures when the 



normal layer is a noncentrosymmetric magnetic metali^. 
Josephson ip junctions demonstrate unusual properties 
and may serve as phase shifters in the superconducting 
(SC) electronics circuitsi^. 

The 7r((/?) junctions described above utilize an intrin- 
sic phase change due to the peculiarities of tunneling 



through the ferromagnetic layer or/and superconduct- 
ing wavefunction symmetry. An alternative approach is 
to produce a phase shift across the junction using flux- 
or current-biasing. Examples are tools based on trap- 
ping fluxoids in a mesoscopic ring incorporated into dc- 
SQUID^^ or Josephson junction with the additional cur- 
rent injector-extractor pair which creates an arbitrary 
discontinuity of the Josephson phase differenc e^^i^^ . Here 
we suggest to use small ferromagnetic particles to cre- 
ate the phase-biased Josephson system. An arbitrary 
phase drop across the junction is shown to be caused by 
a small-scale phase difference inhomogeneity induced by 
the particle. 

Let us briefly remind of the basic mechanisms which 
could provide a strong phase variation along the con- 
tact on a scale which is smaller than the Josephson 
penetration depth. First of all, the natural source of 
a phase inhomogeneity is an Abrikosov vortex (AV) 
pinned in the SC electrodes. Even a single misaligned 
AV, trapped in the SC electrodes perpendicular to the 
junction plane, is known to modify strongly the crit- 
ical current and the current-voltage characteristics of 
Josephson junction a^°i^^'^^i^^'^'* . Next, — tt disconti- 
nuities in the phase difference along the barrier appear 
in a Josephson junction composed of alternating and 
TT mini-junctions^^ (the zigzag junctions between high- 
Tc and conventional superconductors^^ as well as SFS 
junctions with a steplike thickness of the ferromagnetic 
interlayer— ) . Finally, arbitrary phase inhomogeneity can 
be formed with a current injection into Josephson junc- 
tion on a scale smaller than the characteristic Joseph- 
son lengtbiSiii^. The presence of such phase singulari- 
ties results in the unusual CPRf^, an anomalous non- 
Fraunhofer IdH) dependence^, and spontaneous gener- 
ation of fractional Josephson vortices at the boundaries 
between and tt regions?^'^^. 

Another method for creating of a controlled phase 
inhomogeneity in Josephson junctions has been pro- 
posed and successfully realized recently. This method 
is based on the interaction of a Josephson contact with 
the nonuniform magnetic field of submicron ferromag- 
netic particles located close to the junctior^i^i^. These 
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experiments have demonstrated a essential dependence 
of the critical current 1^ on the magnetic state of the 
particles. This means that the transport properties of 
Josephson contacts can be effectively controlled by the 
magnetic field of the small particles. Experimentally 
detected additional maxima in the field dependence of 
the critical current Ic{H) have unambiguously indicated 
commensurability effects between a periodic distribution 
of the Josephson phase difference created by the particles 
and the scale of the phase modulation induced by an ap- 
plied magnetic field H--^^ . While the macroscopic com- 
mensurability effects have already been demonstrated in 
such hybrid fcrromagnet-superconductor (FS) systems, 
the insight into the current-phase relation is still lack- 
ing. 

In this paper we study theoretically the hybrid FS sys- 
tem consisting of a Josephson junction coupled with a 
single magnetic dot as it is shown in FiglT] We demon- 
strate that the phase shift in the ground state of the 
Josephson junction depends on the magnetic state of the 
particle. We discuss the possibility to reahze of Joseph- 
son (p-junction based on such hybrid FS structure. The 
paper is organized as follows. In SecHIl we introduce a 
model used to explain the appearance of a nonuniform 
phase-difference distribution in Josephson junction cou- 
pled with a magnetic dot. In Sec lIIIl we study the ground 
state of this hybrid system which is characterized by the 
finite phase difference drop (po across the junction. In 
Section [I VI we summarize our results. 

II. JOSEPHSON PHASE MODULATION 
INDUCED BY A MAGNETIC PARTICLE 

We consider a generic example of the FS hybrid system 
consisting of a short square {W x W) Josephson junction 
and an elongated magnetic particle on its top electrode. 
The junction is formed by overlapping two long SC strips 
(L > W) of thickness d < A (top) and D >2X (bottom) 
as shown in Figllji. The single-domain magnetic dot with 
in-plane magnetization M is separated from the top SC 
electrode by a thin insulating layer, which prevents the 
proximity effect. The interaction between the junction 
and the particle may be provided by the magnetic field 
generated by the dot and supercurrents. For the sake 
of simplicity, we consider here only the case of a rather 
small junction with 

A<ty<A, Aj, (2) 

where A and Aj are the London penetration depth and 
the Josephson penetration depth, respectively, and A = 
A^/d is the thin-film screening length. The gauge invari- 
ant phase difference across the junction is given by the 
expression 

t 

^{r)^e\r)-e\r) + ^J dz Mr) . (3) 

b 



Here ^''(r) and 9*(r) are distributions of the phase of 
Cooper wave function in the bottom and the top SC elec- 
trodes, respectively, r = {x, y) is a vector in the junction 
plane, Az is the normal to the junction plane component 
of the magnetic vector potential A = A|| + A^zq, and 
<i>o — Trhc/ I e I - is the flux quantum. In a thin SC strip 
with W ^ A, the self-field of the sheet current can be 
disregarded, and the top electrode is assumed to be trans- 
parent to the magnetic field of the particle. In its turn, 
this magnetic field partially penetrates into the bottom 
superconductor and induces in-plane screening Meissner 
currents in it. Since the top film thickness d is assumed 
to be small it is reasonable to apply the gauge A^ — 
to avoid the strong singularity of the vector potential A 
in the limit d ^ 0. For this gauge the phase difference 
across the junction (j>{r) ^ is determined by the distri- 
butions 6** (r) and 6^ (r) of the phase of the wave functions 
only. 

In the absence of vortex lines trapped in the electrodes 
the ground state of the system can be described by the 
uniform phase 0* = 0'' = and screening Meissner cur- 
rents are determined by the the in-plane component of 
the vector potential Ay. So, the gauge invariant phase 
difference ([3]) across the junction equals zero: (j){r) = 0, 
and the magnetic field of the particle does not modify 
the critical current of the Josephson junction. 

Vortex (antivortex) , if appears, must be located near 
the negative (positive) pole of the magnetic dotj^, as 
shown in Fig. la. According to the concept proposed in 
Reffsj. a pair of vortices of opposite directions pierce the 
top electrode of the junction. As before, we put 6^ = Q 
due to the absence of vortex lines in the bottom electrode. 
Suppose that the pair size (i.e., the vortex-antivortex dis- 
tance) |r„ — TqI = 2a is rather large compared with the 
superconducting coherence length ^, then the electrody- 
namic mechanism based on the spatial dependence of the 
gauge invariant phase difference 

Hr) = -0\r) (4) 

is dominant'^^'^^. Here r„ — (xy, y^) and = (xa, ya) 
are the vortex and antivortex positions, respectively. 



A. Basic equations 

As a next step we should find the phase difference 
distribution 0(r) over the junction area. The small 
size of the junction W <^ Xj means that self-field ef- 
fects of the Josephson current can be neglected com- 
pared to the in-plane currents. In this case, the phase- 
difference distribution (f)(r) obeys the two-dimensional 
Laplace equation'^^''^^: 

A(/. = + d^(j) = 0. (5) 

In the presence of trapped Abrikosov vortices the top 
electrode of the Josephson contact becomes a multiply 
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connected domain. So, it is necessary to take into ac- 
count the topological singularities of the phase distribu- 
tion 6'*(r) which are caused by presence of the vortex- 
antivortex pair'^^. According to the Eq. ^ we obtain 
the following condition for 4>{r): 



curl^(V0) = 27r[(5(r 



(6) 



which fixes the circulation around the singularities. 

For W <^ A, the in-plane sheet current density in the 
top electrode 



'87r2A 



27r 



(7) 



is determined mainly by the phase gradient V0* induced 
by the trapped vortices rather than the vector potential 
Aj[. Indeed, the term with A|| in the Eq. ([7]) is of the 
order of 2tt/A < |V6''| 2tt/W and can be neglected. It 
means that the sheet current density ([7]) is determined by 
the gauge invariant phase difference ^(r): cr' VS* ^ 
V(j). At the edges of the top stripe {x = 0,W ; y = 0) 
the normal component of the sheet current at vanishes, 
and the Eq. (O must be supplemented with the following 
boundary conditions: 



dx 



0, 



x=0,W 



0. 



(8) 



Finally, a local phase inhomogeneity due to the presence 
of a vortex-antivortex pair has to vanish at distances 
larger as compared with the pair size 2a. For the sake 
of simplicity the top electrode is assumed to be a semi- 
infinite SC strip, and the condition 



{x, y) = 0, for y 



(9) 



has to be satisfied. Thus, the equations (O,® and the 
boundary conditions (O,® describe the phase difference 
distribution which is induced by a vortex-antivortex pair 
trapped in the top electrode of the contact. 

At a constant value of the critical current density jc 
and the standard sinusoidal form of CPR the ground 
state of this junction in the absence of a supercurrent 
corresponds to a minimum of the Josephson energy 



+ 0(r)) , 



(10) 



where the integral is evaluated over the junction area 
Sj : {0 < x,y < W), and Ic = jcSj. The current-phase 
relation for the junction 

Ij{^)^jcjdrsm{^ + ^{r)), (11) 



reveals the shift which depends on the phase difference 
distribution 0(r). 

By minimizing the Josephson energy pop with respect 
to if one can find the equilibrium distribution 0e(r) = 



(/>(r) -I- ifo of the gauge invariant phase difference for the 
junction containing the vortex-antivortex pair, where the 
additive constant 



ipo = -arctan (S'^/C^) 



(12) 



determines the fixed phase shift between the top and bot- 
tom electrodes away from the junction area. The coeffi- 
cients and in (|12p depend on the distribution of the 



phase difference 
vortices: 



— dr cos 0(r) , S^— I dr sin 0(r) 



(r) induced by the trapped Abrikosov 



(13) 



Varying the magnetic state of the particle one may con- 
trol 0(r) changing the vortex(antivortex) position, and, 
thereby, tune the average phase difference ip^ across the 
junction. 



B. Phase Difference Distribution 

The the gauge invariant phase difference (/)(r) due to 
the presence of a single vortex-antivortex pair trapped in 
the top electrode can be calculated as follows. First of all, 
let's neglect the edge effects ([5]). The distribution 0(r), 
which satisfies the Laplace equation ([5]) and provides for 
the required phase circulation ([6]) around the singularities 
r^^Q, can be written as a superposition of contributions 
from two point opposite vortices: 

(/.(r)=0p(r), c^p{v) = Oaiv) ^ e,{v) . (14) 

The phase 0«,a(r) of the wave function describing a point 
vortex (antivortex) is determined by the polar angle spec- 
ifying the direction from the position of the vortex axis 
Vy^a to the reference point r (see Fig. [TJd): 



6'u,a(r) = arctan 



V - Vv, a 



X 



X,: 



(15) 



It is evident that (j>p{r) for |r| ^ |rt,,a| and, thus, the 
distribution (fT4|) . (fTS]) satisfies the condition ([9]). Figure[2] 
illustrates schematically the distribution of the phase dif- 
ference (j>p{r) created by a pair of opposite vortices. The 
dark central area between the vortices shows the region 
where 7r/2 < (j)p < 3tt/2 and cos0p < 0. This do- 
main provides with an additional positive contribution to 
the Josephson energy (fTU)) and, as a consequence, such 
state of the junction appears to be an energetically un- 
favourable. The energy excess associated with the pair 
presence grows as the intervortex distance |r„ — ra| in- 
creases. 

To take into account the boundary conditions ([8]) the 
solution 0(r) may be written in the following convenient 
form: 



(16) 



4 



Here </'p(r) = 0'^{r) — 0'^{r) is the phase distribution cre- 
ated by the vortice images (see Fig. [Hd) 



in the infinite stripe ( < a; < W, \ y\ < oo ) with the 
following boundary conditions at the stripe edges x = 
and X = W 



0'^ jj(r) = arctan 



y + Vv.a 



(17) 



and '0(r) is the solution of the Laplace equation 

AV = (18) 



x=0,W 



-dx (0p + 



(19) 



■ip{x, y) = 



du 



du 



yv{u^ -pI) 



{pi + u2)2 _ Ay^u^ (p2 -\- u'^Y — ^y^u^ 



In 



In 



cosh TT 



cosh TT 



_ y-u 
w 

y-u 



W 



cos(^-^) 



+ cos(7r — ) 

^ w' 



(20) 



where 

2 _ 2 , 2 2 _ /Tjr _ \2 I 2 

^u,a — ■^v,a ' yv,a t Pv,a — I" ^v.a) + yv,a ■ 

The expressions (I14l)-(I17l),(l20l) determine the distribu- 
tion of the gauge invariant phase difference 0(r) created 
by the pair of opposite point vortices trapped in the thin 
top electrode of the Josephson junction. This phase dis- 
tribution 0(r) is used to calculate the Josephson energy 
and the total current through the contact pT|) for 
different positions of vortex and antivortex. 

III. GROUND STATE OF JOSEPHSON 
JUNCTION COUPLED WITH A MAGNETIC 
PARTICLE 

Now we proceed with the calculations of the ground 
state of the Josephson junction which depends on the 
size and orientation of the vortex-antivortex pair induced 
by the magnetic particle. We restrict ourselves to the 
case of zero homogeneous external magnetic field. For 
simplicity, the vortex and antivortex are assumed to be 
placed symmetrically with respect to the center of the 
junction xq = yo = W/2: 

Xy + Xa ^W, yv +ya^W . 

Figures [3] shows the dependence of the average phase 
difference (/jq ([12) and the Josephson energy Ej{ipo) pU]) 
on the location of the vortices. As an example we con- 
sider the change in the phase fo due to the rotation of 
the vortex-antivortex pair with respect to the midpoint 
xo = yo = W/2. In this case the location of the vortices is 
determined by the intervortex distance 2a and the angle 
of the pair rotation 

. ( yv-ya\ 

a — arctan 1 

\ Xv Xa J 



relative to the direction of x-axis (see Fig. [T]). The range 
of the change in tp^ depends on the pair size 2a. For 
2a <^W (Fig.[3l curve 1) the phase inhomogeneity occu- 
pies a small part of the junction area, and the presence 
of the vortex-antivortex pair does not affect the junc- 
tion properties essentially. The value of (f^ varies weakly 
round about the point lpq — 0, and the junction demon- 
strates mainly the properties of a conventional junction: 
the ground state corresponds to the almost zero phase 
drop across the junction. An increase in the pair size 
leads to forming a strong phase modulation 4'{v), and 
the average phase difference Lp^ can take practically any 
value between — 7r/2 and 7r/2 in dependence on the angle 
a (Fig. [3l curve 2). This is accompanied by the growth 
of the Josephson energy Ej due to the expansion of the 
domain where cos(/)(r) < 0. The further increase in the 
pair size (2a > W/2) leads to the additional tt shift of 
the average phase difference (po (Fig.[3l curve 3), and the 
junction exhibits the switching into the new ground state 
which is specific for tt junctions: the value of (po oscillates 
about the point ^pq — tt. 

Figure [4] shows the dependences of the average phase 
difference ipo and the Josephson energy of the ground 
state Ej° vs the intervortex distance 2a for the fixed ori- 
entation of the vortex-antivortex pair a = tt/2. For a 
small intervortex distance 2a, the domain, where the en- 
ergetically unfavourable phase difference n/2 < 0(r) < 
37r/2 exists, occupies the central part of the junction. 
The size of the domain is small, and this part of the 
junction gives a small additional positive contribution 
to the Josephson energy Ej. As a result, the conven- 
tional ground state realizes, which is described by the 
zero value of the average phase difference tpo. The in- 
crease in the intervortex distance 2a provokes an expan- 
sion of the energetically unfavourable domain, and, con- 
sequently, the rise in the Josephson energy Ej of the 
state corresponding to a choice of (po = 0. Let's in- 
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troduce the complementary phase difference distribution 
(/)* (r) which differs from the initial one ^(r) by the n shift: 
(f>*{r) = (f>{r) + TT. In the absence of trapped vortices the 
complementary state is unrealizable, since it corresponds 
to a maximum of the Josephson energy (jlO[) . The gener- 
ation of the vortex-antivortex pair results in a decrease 
of the Josephson energy Ej of the complementary state 
(f)*{r). The crossing of the curves E'j{2a) and Ej{2a) 
occurs at 2a = 2a* ~ 0.56 VF, so that the switch be- 
tween the conventional state (ipo — 0) and the tt shifted 
one {ipo = T^) takes place. As a result, the ground state 
of the junction changes for sufficiently large intervortex 
distances 2a > 2a*. Certainly, the switching point 2a* 
depends on configuration of the vortex-antivortex pair. 
Such a crossover between and tt states manifests itself 
as a TT shift of the superconducting phase between the 
electrodes of the junction. 

In Fig. Owe present examples of the simulations of the 
equilibrium phase difference 4)(r) + tpo for the junction 
containing the vortex-antivortex pair. For a small inter- 
vortex distance (2a < 2a*), the conventional ground state 
(fiQ — realizes (Fig. [5^). This state corresponds to the 
zero phase drop across the junction in the absence of a su- 
percurrent. With the increase in the intervortex distance 
(2a > 2a*), the additive phase shift (po = tt decreases 
in part the Josephson energy corresponding the central 
domain between the vortices (Fig. [5b). As a result, there 
appears a tt shift in the phase of the superconducting 
order parameter across the junction. 

IV. SUMMARY 

In conclusion, we have studied theoretically the prop- 
erties of the hybrid system consisting of a short Joseph- 
son junction located in a nonuniform field of the mag- 
netic particle. The effect of the particle on the junction 
is shown to be strong due to the formation of the pair 
of oppositely directed Abrikosov vortices which pierce 
the top electrode of the junction. From an experimen- 
tal point of view the vortex-antivortex pair can be cre- 
ated by cooling the junction through transition temper- 
ature Tc in the dipole magnetic field of a magnetic dot'^''. 
As a result, one should use single-domain magnetic dots 
with in-plain magnetization such as elongated submicron 
Co dots. For parameters of magnetic dots taken from 



the experimentSf22, (the saturation magnetization Mg ~ 
800 Oe, lateral dimensions ~ 650 (easy axis) x 250 nm^ 
and a thickness of 50 nm) we can estimate typical 
dimension of the Josephson junction as ^ 1 /im. 
Magnetostatic calculations show that the stray field of 
both poles of the dot creates a (positive or negative) 
fiux $s > $0 through the top SC electrode, and the 
criterion^i^^ for the nucleation of the vortex-antivortex 
pair is satisfied. Since the flux $s rapidly decreases when 
a distance between the magnetic dot and a superconduc- 
tor becomes comparable with the lateral sizes of the dot, 
the Meissner state is more energetically favorable in the 
thick bottom electrode upon cooling through Tc. 

The vortex-antivortex pair trapped in the top elec- 
trode causes the inhomogeneity of the gauge invariant 
phase-difference on scales that are significantly smaller 
than the Josephson length. The pair size and orienta- 
tion are determined by the size and magnetization of the 
particle. We have calculated the corresponding distri- 
bution of the Josephson phase difference over the junc- 
tion and have studied the dependence of the Josephson 
energy Ej on the vortex configuration. It was shown 
that the ground state of the junction corresponds to the 
nonzero drop of the average phase difference across the 
junction ipQ. The equilibrium value of the phase differ- 
ence ipo depends on the size and orientation of the vortex- 
antivortex pair, i.e. on the magnetic state of the particle. 
This demonstrates the possibility to realize the tunable 
current-phase relation 1^ = Ic sin(c/? -I- ipo) by changing 
the magnetization of the particle. Finally, note that such 
hybrid FS system incorporated into a superconducting 
ring should induce spontaneous currents and may serve 
as a natural phase shifter in the superconducting elec- 
tronic circuits. 
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FIG. 1: (Color online) (a) Schematic diagram of the Josephson contact with the ferromagnetic particle on the top electrode. 
The junction area {W x W) occupies only part of the superconducting electrodes {L ^ W). The vortices indicated by vertical 
arrows are located near the opposite poles of a uniformly magnetized particle. The inset schematically shows the structure 
of the stray field of the ferromagnetic particle and the pair of oppositely directed Abrikosov vortices which pierce the top 
superconducting electrode, (b) The picture shows the location of a vortex (r„) and antivortex (va) trapped in the top electrode 
of the contact. The phase difference (pp is specified by the angle between the directions from the reference point r to the points 
of the vortices location rv,a- The thin dashed lines shows vortex(antivortex) images, which provides zero normal currents at 
the edge y — 0. 
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FIG. 3: (Color online) Dependence of the average phase difference <^o (a) and the Josephson energy i?™ of the ground state; 
(b) on the angle of the vort ex-ant i vortex pair rotation a for different values of the pair size 2a: \. 2a = Q.\W\ 2. 2a = 0.414^; 
3. 2a = 0.8W (So = hlc/2e). 



0.0 0.2 0.4 0.6 0.8 1.0 



2a /W 



FIG. 4: Dependence of the average phase difference tpo and the Josephson energy E'j° of the ground state on the size 2a of 
the vortex-antivortex pair (a = 7r/2). The dashed hnes show the dependence E^{2a) and Ej{2a) for ipo — and cpo = tt, 
respectively; Ej = min{£;5, EJ} (Eo = ft/c/2e). 





FIG. 5: Distribution of cosine of the phase ip — <j!>(r) + ifo for different values the pair size (intervortex distance) 2a: (a) 
2a — 0.5 W ; (b) 2a = 0.6 H^. The edge of the junction y — W is marked by the dash line. 



